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Abstract 

A semi-microscopic self-consistent quantum approach developed recently to describe 



, the inner crust structure of neutron stars within the Wigner-Seitz (WS) method with 

CN I the explicit inclusion of neutron and proton pairing correlations is further developed. In 

' this approach, the generalized energy functional is used which contains the anomalous 

. term describing the pairing. It is constructed by matching the realistic phenomenological 

functional by Fayans et al. for describing the nuclear-type cluster in the center of the WS 
cell with the one calculated microscopically for neutron matter. Previously the anomalous 
part of the latter was calculated within the BCS approximation. In this work corrections 
I to the BCS theory which are known from the many-body theory of pairing in neutron 

matter are included into the energy functional in an approximate way. These modifications 
■ have a sizable influence on the equilibrium conflguration of the inner crust, i.e. on the 

proton charge Z and the radius Rc of the WS cell. The effects are quite significant in the 
' region where the neutron pairing gap is larger. 

\o 

O ' PACS. 26.60.-I-C Nuclear matter aspects of neutron stars - 97.60. Jd Neutron stars - 21.65.-|-f 
Nuclear matter - 21.60.-n Nuclear structure models and methods - 21.30.Fe Forces in hadronic 
systems and effective interactions 

O 

F; ; 1 Introduction 

> . 
•1—1 . 

^ ■ Recently, we have developed a semi-microscopic self-consistent approach to describe the neutron 
star inner crust with the explicit inclusion of neutron and proton pairing correlations PJ E] . 
The inner crust is the part of the neutron star shell with subnuclear densities, O.OOlpo ^ 
p < O.Spo, where Pq is the normal nuclear density. It is a crystal system consisting mainly 
of spherically symmetrical nuclear-like clusters immersed in a sea of neutrons and virtually 
uniform sea of electrons. The quantum self-consistent description of the inner crust goes back 
to the classical paper by Negele and Vautherin (N&V) |3] who used a kind of energy functional 
method combined with the Wigner-Seitz (WS) method to describe the crystal structure effects 
in an approximate way. Within this approach, for a fixed average nuclear density p, the energy 
functional of the system is minimized for the spherical WS cell of the radius Rc- A cell contains 
A = {4:7i/3)R^p nucleons, specifically Z protons and N = A — Z neutrons. As far as the system 
is electro-neutral, the number of electrons per a cell is equal to Z. For a mature neutron 
star which can be considered at zero temperature and neutrino free, the /5-stability condition 
consists in equality of the neutron chemical potential /i„ to the sum of the proton and electron 
ones, Pp+pe- For a wide region of p, the minimization procedure was carried out in j3jj for 
different values of Z and Rc resulting in the equilibrium configuration (Z, Rc) for the density 
under consideration. 
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It should be noted that the pairing effects were not taken into account in |3j. Thus, the 
approach developed in |2j could be considered as a generalization of the N&V method with 
allowance for pairing effects. Although the contribution of the pairing to the total binding 
energy of the system under consideration is rather small, it turned out that it may to change 
the equilibrium configuration {Z, R^) significantly. 

To involve the pairing effects in a self-consistent way, we use the generalized energy func- 
tional method |3] which incorporates in a natural way the pairing into the original Kohn-Sham 
jS] method. In this approach, the interaction part of the generalized energy functional is the 
sum of the normal component, £^norm(p), and the anomalous one, ^an(p, J^)- They depend on the 
normal densities {pn, Pp) and the anomalous ones, {i>n, i^p)- Just as in the Kohn-Sham method, 
the prescription m* = m holds to be true. The so-called Superfluid LDA method suggested 
recently is rather close to the method by S. Fayans et al. The main difference between the 
two approaches is in the form of the density dependence of the anomalous term of the energy 
functional. 

In P3I2] the semi-microscopic energy functional was constructed with matching at the nuclear 
cluster surface the phenomenological nuclear functional E"^^ by Fayans et al. inside the cluster 
to a microscopic one, £^™, for the neutron environment. The normal part of was found in [7] 
within the Brueckner approach with the Argonne vig potential [S] . The anomalous component 
of was calculated in |2j within the Bardeen-Cooper-Schrieffer (BCS) approximation for 
neutron matter, again with the use of the vig potential. 

It is well known that the BCS approximation overestimates the gap value A„ in neutron 
matter. Various many-body corrections suppress the value of Aj^*-'^ significantly. Although up 
to now there is no consistent many-body theory of pairing in neutron matter, there exists a 
conventional point of view (see e.g. P) that the BCS gap value is suppressed, 

An(A;, k^) = fm-b{k, fcF)Af kp), (1) 

by a factor /m_b(/c, /cf) which is between 1/2 and 1/3. The only exception seems to be the 
work in ref . ^U] . In this paper, we use a simple model for the many-body corrections in which 
the suppression factor fm-h{k, kp) is supposed to be momentum and density independent. This 
ansatz is essentially similar to that used in |llj for considering the structure of a superfluid 
vortex in neutron matter. We modify the anomalous component of the energy functional of 
d 12 by introducing the constant suppression factor. We examine two versions of this model, 
the P2 model (/m-b = 1/2) and the P3 one (/m-b = 1/3)- In this notation, it is natural to 
name the BCS approximation (/m-b = 1) as the PI model. We expect that the real truth is 
somewhere between the P2 and P3 models. 

One more remark should be made before going to the body of the article. As it was 
found recently ^2); there are internal uncertainties inherent to the WS method applied to 
the neutron star inner crust. They originate from the kind of the boundary conditions for 
the single-particle functions used in the WS method. There are two kinds of such boundary 
conditions which a priori seem equivalent. As it turned out, the corresponding predictions 
for the equilibrium configurations {Z, Rc) are in general different. As a rule, the difference is 
not large, corresponding to variation of Z by 2 - 6 units and of Rc, by 1 - 2 fm. However, 
sometimes strong changes in the neutron single-particle spectrum arise which influence the 
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solution of the gap equation significantly. Therefore, for each model under consideration, we 
carried out calculations for both kinds of boundary conditions. 

2 Modification of the BCS anomalous part of the Gen- 
eralized Energy Functional 

The ansatz of IH |2] for the complete energy functional consists in a smooth matching of the 
phenomenological and the microscopic functionals at the cluster surface: 

^(p.(r), Urir)) = S^\p,{r), Mr)) + {S^\pr{r), M^)) - ^"'(Pr(r), z/.(r))) (1 - F^r)) , (2) 

where t = n,p is the isotopic index and the matching function F^ir) is a two-parameter Fermi 
function. The latter is taken to be the same for the normal and the anomalous components of the 
energy functional, with the diffuseness parameter dm=0.3 fm and the matching radius Rm which 
should be chosen anew in any new case, in accordance with the equality of Pp{Rm) = O.lpp(O). 
For such a choice, practically all the protons are located inside the radius Rm- Therefore, 
the matching procedure concerns,in fact, only neutrons, protons being described with the pure 
phenomenological nuclear energy functional. In practice, we use in Eq. Q an approximation in 
which only neutron components of the microscopic and phenomenological functionals are taken 
into account in the second term containing the difference of (£™ — S^^). 

Following to m 12] ; we use for the microscopic part of the normal component of the total 
energy functional Q the one calculated in [7] for neutron matter with the Argonne Vig potential. 
Its explicit form could be found in the cited articles. Here we concentrate on the anomalous 
part of the energy functional which will be modified in comparison with that of 

The anomalous part of the energy functional used in ^1 |2j has the form: 

^.n = lT.Krir)Wr{r)\\ (3) 

where Vf^^ is the density dependent effective pairing interaction. 

The matching relation (j2)) for the anomalous part of the energy functional leads to the 
analogous relation for the effective pairing interaction: 

Vf (r) = V!^{p{r))FUr) + VMr)){l ' FUr))- (4) 

The isotopic index r is for brevity omitted. 

We shall use the same phenomenological effective pairing interaction Vg^ as in I2] , there- 
fore we omit here its explicit form. Let us note only that it has a density dependent coordinate 
delta-function form of |4.. The explicit form of the density dependence |3j was modified a 
little in pQ in accordance with (@)). As to the microscopic effective pairing interaction, it was 
calculated in PP microscopically within the BCS approximation with the same Argonne force 
Vi8 as the normal part of the energy functional . In this paper, we modify this procedure to 
take into account approximately the many-body corrections to the BCS approximation. Let us 
first repeat the BCS procedure. 
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The microscopic part of the effective pairing interaction, V^{r), should be found for the 
model space 5*0 under consideration which is limited with the energy Eq for the single-particle 
spectrum. For a fixed value of the neutron density Pn{r), it is defined via the gap equation in 
homogeneous neutron matter with the density p = Pn{f)- We start from the BCS approximation 
in which the gap A is expressed directly in terms of the bare NN potential f (fc, k') in the 15*0 
channel: 

m = -E^K*^.*')^. (5) 



where E{k) = y (e^ — + ^^(A;), Sk = k'^/2m + Un, Un is the value of the neutron matter 
potential well. In terms of the effective pairing interaction, the gap equation looks analogously, 
but the integration in the momentum space is limited within the model space So'- 



where ko = ^^2m{Eo + /i„ - Un)- 

In the BCS approximation the relation between the effective pairing interaction and the 
bare NN potential is obvious: 

V^iik, k') = vik, k') - <k, ^i):TT^Vrg\fci, k')- (7) 

The effective pairing interaction entering Eq.(jni) depends explicitly on momenta, which 
corresponds to a non-local force in the coordinate space. In view of very simple local form of 
the phenomenological effective pairing interaction V^g in Eq. for matching it is necessary 
to simplify the microscopic partner V^' to a local form, too. The simplest way is, for a fixed 
value of p under consideration, to define it as a fc-independent average value of the effective 
pairing interaction in Eq. © which yields the same value A(A;f) as the exact effective pairing 
interaction: 



where kp is the local Fermi momentum kp = (Svr^p)^/^, and EQ{k) = \J{ek — PnY + A2(A;p). 
For the Fermi momentum k-p under consideration, the microscopically calculated value of A(A;f) 
is the only input of Eq. (jH)) to find the effective pairing interaction V^{kp)- In the BCS (or PI) 
model, one uses A^^^(/cf) in Eq. (jH)). In the P2 or P3 models the value of ^{kp) is found from 
Eq. (d} with /m-b=l/2 or 1/3, correspondingly. Then it should be substituted into Eq. (jH}. 
The resulting values of the effective pairing interaction in neutron matter for the case of the 
PI (BCS), P2 and P3 models are displayed in Fig. 1. 

Let us now discuss the problem of the boundary conditions in the WS method mentioned 
above. For the case of the BCS approximation, it was examined in ^21- Application of the vari- 
ational principle to the energy functional under consideration for a WS cell results in the set of 
the Shrodinger-type equations for the single particle neutron functions (pxij) = i?„;j(r)$/jm(n). 
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Figure 1: The effective pairing interaction in neutron matter for tlie PI , P2 and P3 models 
(solid, dashed and dotted line, correspondingly). 
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with the standard notation. The radial functions Rnijir) obey the boundary conditions at the 
point r = Re- There exist different kinds of the boundary conditions. N&V used the following 
one: 

Kz,(r = = (9) 

for odd /, and 

K,^.(r = R,) = 0, (10) 

for even ones. In it was denoted as BCl. An alternative kind of the boundary conditions 
(BC2) was considered also there, when Eq. Q is valid for even / whereas Eq. ()10j). for odd ones. 
As the analysis of |12| for the BCS case (PI model) has shown, some predictions of the two 
versions of the boundary conditions (BCl versus BC2) are in general different. In this paper, 
we will perform the similar analysis for the P2 and P3 models. 

3 A brief summary of the pairing effects in the case of 
the BCS approximation 

Calculations of jH |2| were carried with the use of the BCS approximation for the neutron 
matter pairing and the N&V version of the boundary conditions (BCl, in our notation). As it 
turned out, the pairing correlations influence the equilibrium values of {Z, R^) significantly. In 
the paper [T2], for the case of the BCS approximation (i.e., the PI model), the dependence of 
the ground state properties of the inner crust on the kind of the WS boundary conditions was 
examined. The calculations were carried out directly with the two kinds of boundary conditions, 
BCl and BC2, and the results were compared with each other. To make the analogous analysis 
for the P2 and P3 models in the next sections more transparent, we cite here some results of 
jl2j . In particular. Fig. 2 shows the values of the binding energy per a nucleon, E-q, for the PI 
model calculated with the two versions of the boundary conditions. 

We see that in the case of very small density, /cf=0.2 fm^^, which is nearby the neutron 
drip point, the predictions of the BCl and BC2 versions are practically identical. At increasing 
density, with fcp > 0.6 fm~^, the uncertainty in the equilibrium value of Z is between 2 and 
6 units, with the largest values at the largest k-p. The uncertainty in the value of Rc is, as a 
rule, about 1 fm and only for k-p=l.l fm~^ it turns out to be about 2 fm. Referring to [21 
for more details, we present in Table 1 the main ground state characteristics of the inner crust. 
There are two lines for each value of k-p. The first one is given for the Z value corresponding 
to the minimum of the binding energy E-q in the BCl case, the second one, for BC2. The 
only exception is /cf=0.2 fm~^ where these two values of Z coincide. It should be stressed 
that, as a rule, the value of these uncertainties is smaller than the variations of the equilibrium 
configuration (Z, Rc) connected with the pairing effects [2]. One more point which is important 
for the analogous consideration within the P2 and P3 models is as follows. For all the values 
of kp which were investigated, the relative position of the local minima of the functions Ep,{Z) 
is always similar for the BCl case and the BC2 one. As the result, the corresponding absolute 
minima are rather close to each other. 

Table 2 collects some important characteristics of the gap for all the density values under 
consideration. The relative position of two lines for the same value of kp is the same as 
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Figure 2: Binding energy per a nucleon for various in the BCl case (solid circles connected 
with the solid lines) and the BC2 one (open circles connected with the dotted lines). 



7 



Table 1: Comparison of characteristics of equilibrium configurations of the WS cell for two 
different kinds of the boundary conditions in the case of the PI model (BCS approximation). 



kp, 
fm^^ 


Z 


Rc, fm 


Eb, MeV 


Atn, MeV 


BCl 


BC2 


BCl 


BC2 


BCl 


BC2 


0.2 


52 


57.18 


57.10 


-0.9501 


-0.9483 


0.1928 


0.1942 


0.6 


58 
56 


37.51 
36.97 


37.48 
36.95 


2.1516 
2.1563 


2.1596 
2.1572 


3.2074 
3.2173 


3.2226 
3.2193 


0.7 


52 

48 


32.02 
31.16 


32.04 
31.14 


2.7908 
2.7924 


2.7989 
2.7856 


3.9876 
4.0069 


4.0107 
3.9873 


0.8 


42 
11 


26.90 
27.29 


26.91 

27.30 


3.4373 
3.1135 


3.4471 

3.1319 


4.8454 

1.8553 


4.8561 

1.8198 


0.9 


24 
22 


20.26 
19.87 


20.30 
19.70 


4.1123 
4.1141 


4.1169 
4.1104 


5.7340 
5.7861 


5.7986 
5.7170 


1.0 


20 
24 


16.69 
18.29 


16.90 
18.22 


4.8210 
4.8479 


4.8522 
4.8231 


6.8525 
6.8446 


6.7424 
6.8920 


1.1 


20 
26 


14.99 
16.75 


15.33 
17.08 


5.5765 
5.6677 


5.6733 
5.6100 


7.4288 
7.9680 


8.0446 
8.5398 


1.2 


20 
26 


13.68 
15.21 


13.95 
14.89 


6.4225 
6.6639 


6.6762 
6.4587 


8.5814 
9.0825 


9.1898 
9.3413 
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in Table 1, i.e. the first line corresponds to Z from the equilibrium {Z,Rc) configuration 
for the BCl version. The following notation is used. The asymptotic value of the Fermi 
momentum kp corresponds to the asymptotic value of the density p(r) averaged over the 
interval Rc—b < r < Rc, 6=2 fm. The asymptotic gap value Aas is found as the average of A(r) 
over the same interval. The central gap value A(0) is calculated as the average of A(r) over the 
interval < r < 3 fm. The Fermi average gap Ap is defined as the average value of the diagonal 
matrix element of the neutron gap at the Fermi surface. The averaging procedure involves 10 
levels above and 10 levels below. At last, Ainf means the infinite neutron matter gap value 
found within the BCS approximation for the density p corresponding to the Fermi momentum 
kp, and A^^^ is the same for the Fermi momentum kp. Let us remind that the latter corresponds 
to the average nucleon density under consideration. Obviously, the inequality kp < kp takes 
place because the WS cell contains a nuclear-like cluster in the center with the density which 
exceeds the average one. The difference is especially large in the case of kp=0.2 fm which is 
nearby the neutron drip point. Indeed, in this case almost all the matter is concentrated in the 
central blob. So big difference between the values of Ajnf and A?^f for kp=0.2 fm is explained, 
first, by the big difference between two values of the Fermi momentum and, second, with the 
sharp dependence of the A„ in neutron matter on kp at small kp. 

It is worth to mention that the difference between the asymptotic Aas value and the infinite 
neutron matter prediction Ajnf is a measure of validity of the LDA for the gap calculation 
outside the central nuclear cluster. One can see that, as a rule, the LDA works within 10% 
accuracy, but sometimes the difference is greater which is an evidence of the so-called proximity 
effect. The Fermi average value Ap is usually very close to Aas value. It is explained with the 
fact that the region out of the nuclear cluster, in which the function A(r) is almost a constant, 
contributes mainly to the matrix elements of A nearby the Fermi surface. 

For the case of small and intermediate densities, kp < 1 fm~^, the influence of a particular 
choice of the boundary conditions, BCl or BC2, to the value of Ap or Aas is not essential. As 
the result, the uncertainty in predictions for the gap function A(r) caused by this choice of 
the boundary conditions is also rather small. An example for kp=0.8 fm is given in Fig. 3. 
The difference between any couple of these curves is less than the accuracy of the approach, 
and any of them could be used as a prediction for A(r). To be definite, let us consider the 
"self-consistent" gap function for the BCl version of the boundary conditions as the prediction 
of the WS method for A(r) in the case of small and intermediate densities, kp < 1.0 fm~^. In 
the case of fcF=0.8 fm~^ under consideration it corresponds to Z=42. Such a choice is similar 
to that used in [2 . 

Fig. 4 collects predictions for A(r) in the BCS approximation for these values of kp. In 
accordance with the above agreement, the BCl kind of the boundary conditions is used. The 
value of kp=l fm~^ is included as an optional one as far as in this case the uncertainty is not 
negligible (about 20%) but it is not so big as the one for higher values of kp. 

On the contrary, as it is seen in Table 2, in the case of high densities, kp ^ 1 fm~^, the 
uncertainty in the value of the neutron gap is rather large. As it was shown in ^2], such big 
variations (BCl versus BC2) of the gap value in the WS approximation appear due to the shell 
effect in the neutron single-particle spectrum which is rather pronounced in the case of high kp 
and, correspondingly, small Rc values. We consider this effect as an artifact of the WS method 
which should disappear in a more consistent approach. In T^, we suggested an approximate 
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Table 2: Average gap characteristics in the PI model (BCS approximation). 



fm^^ 


Z 


kp, fm ^ 


A(0), MeV 


Aas, MeV 


Af, MeV 


Ainf, MeV 


MeV 








RP? 


H) 1 


Rr'9 


Rn 


Rr'9 


Rn 


Rr9 


0.2 


52 


0.1156 


0.1095 


0.088 


0.132 


0.042 


0.046 


0.043 


0.058 


0.126 


0.106 


0.40 


0.6 


58 
56 


0.5786 
0.5783 


0.5783 
0.5786 


"1 AHA 

1.464 
1.456 


1.471 
1.428 


1.947 
1.899 


1.899 
1.912 


1.919 
1.893 


1.893 
1.891 


2.321 
2.319 


2.320 
2.321 


2.42 


0.7 


52 
48 


0.6758 
0.6763 


0.6753 
0.6763 


1.665 
1.679 


1.650 
1.648 


2.358 
2.312 


2.288 
2.368 


2.300 
2.290 


2.247 
2.325 


2.680 
2.682 


2.678 
2.682 


2.76 


0.8 


42 
44 


0.7729 


n 779zL 

0.7727 


1 7fi7 

1.747 


1.834 


2.580 


2.679 


2.525 


2.560 


z.ooo 
2.883 


9 889 

2.883 


2.93 


0.9 


24 
22 


0.8694 
0.8725 


0.8693 
0.8664 


1.862 
1.936 


1.664 
1.654 


2.777 
2.677 


2.625 
2.680 


2.636 
2.617 


2.506 
2.544 


2.919 
2.918 


2.919 
2.919 


2.92 


1.0 


20 
24 


0.9499 
0.9612' 


0.9613 
0.9574 


1.249 
1.894 


1.966 
1.504 


2.199 
2.705 


2.635 
2.507 


2.023 
2.519 


2.517 
2.288 


2.800 
2.774 


2.773 
2.782 


2.68 


1.1 


20 
26 


1.0315 
1.0434 


1.0531 
1.0649 


0.996 
1.927 


1.889 
1.296 


1.477 
2.469 


2.411 
2.242 


1.318 
2.280 


2.317 
2.020 


2.550 
2.500 


2.458 
2.408 


2.26 


1.2 


20 
26 


1.1243 
1.1278 


1.1321 
1.1160 


1.556 
0.760 


0.992 
0.991 


1.340 
1.549 


2.017 
0.963 


1.210 
1.249 


1.558 
0.862 


2.113 
2.092 


2.066 
2.163 


1.66 



recipe to avoid this uncertainty for the gap function A(r). This topic will be discussed in more 
detail in the next two sections where the P2 and P3 models are considered with the many-body 
corrections to the BCS approximation taken into account. 

4 Pairing effects in the case of the P2 model 

Let us go to the P2 model in which the many-body corrections to the BCS approximation are 
taken into account in an approximate way outlined in Section 2, with the factor f^_\y=l/2 in 
Eq. for the gap A{k-p) which is the input for finding the effective pairing interaction from 
Eq. dHJ. The calculation scheme itself and the presentation of results are mainly similar to those 
for the PI model in the previous section. In particular, again two kinds of boundary conditions, 
BCl and BC2, are used, and the dependence of the results on the choice is analyzed. Fig. 5 
shows the values of the binding energy per nucleon for the P2 model, to be compared to Fig. 2, 
where the PI model is considered. It is not so detailed as Fig. 2 because, according to the 
experience within the PI model, we limit the analysis to the vicinity of the absolute minimum 
of the function Eb{Z) as found for the PI model. Let us first discuss the results with the BCl 
boundary conditions. In this case, a detailed analysis was made for A;f=0.9 fm~^, kp=l.l fm~^ 
and /cf=1-2 fm~^. Comparison with Fig. 2 shows that, at a fixed value of kp, the positions of the 
local minima in the P2 model are close to those in the PI model. Besides, the relative positions 
of the absolute minimum and the local ones are the same for the two models. Therefore for 
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Figure 3: The neutron gap for A;f=0.8 fm~^, and Z=44, in the BCl case (sohd hues) 

and in the BC2 one (dashed hues). 

other values of /cp we examined only the vicinity of the corresponding absolute minimum in 
the PI model. The only exception is the "suspicious" case of /cf=0.7 fm~^ for which, in the PI 
model, the values of Eb{Z) in two local minima are rather close. The systematic calculation 
(with the step SZ^A, instead of the regular one, 5Z—2) showed that in this case the general 
rule formulated above is also valid. 

Let us consider now the results in the case of the BC2 kind of boundary conditions. Again 
the situation is quite similar to that in the PI model, as confirmed by the detailed comparison 
of the two sets of calculations for the BCl and BC2 versions in the case of A;f=0.9 fm~^. 
Therefore, just as in the PI model, in the case of the BC2 boundary conditions one can limit 
the analysis in the vicinity of the absolute minimum for the BCl one. 

Table 3 presents the main ground state characteristics of the neutron star inner crust within 
the P2 model. It is organized similarly to Table 1. One can see that for all the values of kp under 
consideration, with the exception of kp—O.S fm~^, the equilibrium values of Z in the BCl and 
BC2 cases arc different. Let us now compare the equilibrium Z values in the P2 model with those 
in the PI model (Table 1) for the BCl case. One can see that at high densities, kp > 1.0 fm~^, 
they coincide, being equal to Z=2Q. The maximal difference, dZ=Z^^^^—Z^^'^^=6, occurs in the 
case of A;f=0.7 fm~^. It can be explained with the very fiat dependence of on Z in this case. 
Therefore any change of the calculation parameters, namely, of the gap value in neutron matter 
A„ in the P2 model versus the PI one, can shift the position of the minimum significantly. The 
difference of the WS cell radius values, 5Rc—R''f^^—R^^\ at a given kp appears mainly due to 
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Figure 4: The gap A(r) for A;f=0.6 -j- 1.0 fm ^, the PI (BCS) model, and the BCl version of 
boundary condition. 
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Table 3: Comparison of characteristics of equilibrium configurations of the WS cell for two 
different kinds of the boundary conditions in the case of the P2 model. 



fm~^ 


Z 


fm 


^B, MeV 


/in, MeV 


BCl 


BC2 


BCl 


BC2 


BCl 


BC2 


0.6 


56 
54 


36.85 
36.02 


36.88 
36.04 


2.2837 
2.2838 


2.2842 
2.2839 


3.4755 
3.4660 


3.4818 
3.4629 


0.7 


46 
48 


30.31 
30.98 


30.28 
30.88 


2.9320 
2.9332 


2.9354 
2.9348 


4.2892 
4.3004 


4.2451 
4.2533 


0.8 


40 


25.97 


26.19 


3.5781 


3.5807 


5.0762 


5.0987 


0.9 


20 
26 


18.34 
20.38 


18.60 
20.93 


4.2173 
4.2416 


4.2452 
4.2273 


5.7027 
5.8894 


5.8477 
6.0839 


1.0 


20 
24 


16.56 
18.25 


16.94 
17.80 


4.8641 
4.9411 


4.9427 
4.9013 


7.5247 
7.0171 


6.8272 
6.3298 


1.1 


20 
24 


15.05 
16.56 


15.42 
16.16 


5.5734 
5.7339 


5.7363 
5.6387 


9.0601 
8.2971 


7.9905 
6.9493 


1.2 


20 
26 


13.73 
15.30 


14.90 


6.4175 
6.9244 


6.4566 


8.1824 
8.9542 


10.4879 



the difference in bZ. It is usually of the order of 1 fm, and only in the same case of A;f=0.7 fm~^ 
8Rc exceeds 2 fm. Let us now analyze the influence of the kind of boundary conditions on the 
ground state characteristics in the case of the P2 model. In general, this effect is of the same 
magnitude as in the PI model. Again the uncertainty in the equilibrium value of Z is between 
2 and 6 units and in the value of i?c, about 1 — 2 fm. A more detailed comparison with Table 1 
shows that, as a rule, the effect of the boundary conditions in the P2 model is less than that in 
the PI one, but only a little. For example, at A;f=0.8 fm~^ the equilibrium Z values for the BCl 
and BC2 cases are now equal to each other. However, there is the case of fcF=0.9 fm~^ in which 
the effect under discussion is stronger in the P2 model. It should be noted that in the case of 
fcF=1.2 fm~^ the first line of the table {Z=20) contains empty positions corresponding to the 
BC2 version. This means that in the case under consideration (the P2 model, the BC2 type of 
boundary conditions and Z=20) the WS solution of the type we consider, i.e. the WS cell with 
a nuclear type cluster in the center, is not stable. This is a signal of proximity to the point 
of instability for the phase transition to the homogeneous state. In fact, for such high density 
values corresponding to kp=l.2 fm~^ [TH US] and, maybe, to /cf=1-1 fm~^ [T3] the so-called 
"spaguetti" phase should appear which can not be described within the WS method with the 
spherical symmetry assumed. Therefore our consideration of these kp should be considered as 
optional, and the corresponding results are reported mainly for methodological reason. 

Let us now turn to the analysis of the neutron gap. The main gap characteristics in the P2 
model are collected in Table 4 which is analogous to Table 2 in the previous section. Comparison 
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Table 4: Average gap characteristics in the P2 model 





Z 


kp, im"^ 


A(0), MeV 


Aa., MeV 


Af, MeV 


Ainf, MeV 


MeV 


BCl 


BC2 


BCl 


BC2 


BCl 


BC2 


BCl 


BC2 


BCl 


BC2 


n fi 


56 
54 


0.5797 
n 5783 


0.5790 
n 5781 

yj.o 1 ox 


0.971 
n 980 


0.962 
946 


1.058 
1 063 


1.056 
1 049 


1.059 

1 050 

X .uuu 


1.051 
1 048 


1.163 
1 160 

X . xuu 


1.161 
1 1 59 

X . XUC/ 


1 91 


n 7 


46 
48 


0.6760 
0.6784 


0.6760 
0.6750 


1.057 
1.141 


1.104 
1.044 


1.241 
1.248 


1.244 
1.202 


1.211 
1.235 


1.244 
1.213 


1.340 
1.345 


1.341 
1.339 


1 .oo 


0.8 


40 


0.7691 


0.7752 


1.025 


1.231 


1.264 


1.360 


1.256 


1.346 


1.438 


1.443 


1.46 


0.9 


20 
26 


0.8514 
0.8670 


0.8769 
0.8820 


0.608 
1.113 


1.271 
1.325 


0.834 
1.305 


1.374 
1.323 


0.816 
1.320 


1.401 
1.320 


1.460 
1.459 


1.459 
1.459 


1.46 


1.0 


20 
24 


0.9399 
0.9630 


0.9675 
0.9344 


0.024 
1.465 


1.523 
0.447 


0.041 
1.339 


1.329 
0.644 


0.037 
1.426 


1.453 
0.621 


1.411 
1.385 


1.380 
1.418 


1.34 


1.1 


20 
24 


1.0315 
1.0549 


1.0580 
1.0258 


0.015 
1.439 


1.576 
0.804 


0.026 
1.253 


1.245 
0.649 


0.022 
1.393 


1.424 
0.658 


1.275 
1.226 


1.219 
1.287 


1.13 


1.2 


20 
26 


1.1253 
1.1306 


1.1146 


1.229 
0.177 


0.050 


0.461 
0.305 


0.068 


0.539 
0.309 


0.061 


1.053 
1.037 


1.086 


0.83 



of these two tables shows that the main effect of changing from the PI to the P2 model is quite 
trivial. It is a general suppression of the gap characteristics approximately by a factor two in 
comparison with the PI model (i.e. within the BCS approximation). The ratio A^^'^^ / A^^^^ is 
equal to 1/2 exactly for the quantity A?^^ and approximately for Ainf. For other characteristics, 
the deviation of the ratio from 1/2 is usually a little greater than for Ainf and only in few cases 
it is significant. The Fermi average value Ap is the most important neutron gap characteristic. 
Let us compare its values for the BCl and BC2 versions of the boundary conditions. Just 
as in the PI model, the difference Ap*-^^— Ap*-^^ is very small for A;f=0.6 -^ 0.8 fm~^. For 
kp > 0.8 fm~^ the influence of the boundary conditions on the Ap value becomes rather strong, 
significantly stronger than in the PI model. Especially strong effect appears for kp—l.O im~^ 
and kp=l.l fm~^ where Af almost vanishes in the BCl case, being rather big in the BC2 one. 
At kF=1.2 fm~^, the contrary situation takes place, i.e. Ap vanishes in the BC2 case. 

For the case of kF=l.l fm~^, the A(r) function is displayed in Fig. 6 for the two kinds of 
boundary conditions. The equilibrium value is Z=20 in the BCl case and Z=24 in the BC2 
one. As one can see, the difference between predictions of the two kinds of boundary conditions 
is drastic. The most strong variation of the gap occurs in the case of Z=20. To understand 
the reason of this effect, it is instructive to examine the neutron single particle spectrum ^a- 
It is displayed in Fig. 7 for the BCl version and the BC2 one. The positions of the chemical 
potential /i„ are shown with dots. The two spectra are essentially different. The reason is the 
shift Sex of each A-level going from BCl to BC2 version. The value of the shift is approximately 
equal to one half of the distance between two neighboring levels with the same (/, j), the sign 
of the shift being opposite for even and odd /. The absolute value of the shift is proportional 
to 1/Rl and grows at increasing values of kp. These shifts are shown for two states, 2ji3/2 and 
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Figure 5: Binding energy per a nucleon for various k-p within the P2 model in the case of the 
BCl (sohd circles) and the BC2 (open circles) kinds of the boundary conditions. 
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1^23/2, which are the neighbors of /i„ in the BCl case. One can see that in both cases there is 
a shell type structure with rather wide intervals between some neighboring levels. In the BCl 
case, we deal with a big inter-level space just at Fermi surface, fin being inside. The width of 
this interval exceeds the value of 2Ainf — 2.5 MeV which is characteristic for the gap equation. 
That is why the gap equation in the BCl case has practically zero solution (Ap << Ajnf). In the 
BC2 case, big intervals are situated far from the Fermi surface and do not influence significantly 
the gap equation. Therefore we have a normal situation in this case with Ap — Ai^f . For Z=24 
which is the equilibrium value for the BC2 case, the difference between predictions of the BCl 
and BC2 versions is not so dramatic but also exists. 

Let us return to the PI model (Table 2). One can see that in the case of /cf=1-1 fm~^ and 
Z=20 the Fermi average gap value Ap is suppressed in comparison with the normal one, but 
is not zero. For the solution of the gap equation one can use the same spectrum ex displayed 
in Fig. 7. Indeed, the only difference between PI and P2 models at a fixed value of Z is the 
value of the gap in neutron matter, A^j^j, which is a parameter of the model. In the PI model 
it is about two times larger than that in the P2 model. However, the direct influence of the 
gap on the mean field potential and the single-particle spectrum is negligible. But now the 
value 2Ainf ~ 5 MeV is of the order of the energy interval under discussion. Therefore the 
suppression effect in the gap equation is less than the one in the P2 model. 

The "Shell effect" in the neutron single-particle spectrum was discussed in detail in ^21- It 
was interpreted there as an artifact of the WS approach which does not take into account the 
periodicity of the crystal. It should disappear in a more consistent approach to the neutron star 
inner crust structure with periodical boundary conditions. A recipe was suggested in ^2] for 
improving this drawback and finding the gap A(r) in such anomalous cases with big difference 
of the gap values in the BCl and BC2 versions, one of them being strongly suppressed. It is 
based on the smooth dependence of A(r) on Z in a regular situation and consists in the use 
of A(r) for a neighboring Z with a regular single-particle spectrum at the Fermi surface. In 
the case under consideration. Fig. 6, the solid line (BCl) for Z=2A or the dashed one (BC2) 
for Z=20 correspond to the "normal" situation. The difference between these two curves is 
not greater than 10%, and, within such accuracy, any of them could be used as the prediction 
for A(r) in the case of A;p=l.l fm~^ within the P2 model. As it was discussed in ^21, such 
a recipe is not self-consistent within the WS method but it seems to be reasonable from the 
physical point of view. The arguments were given in this article in favor of the conclusion that 
the solution of the gap equation in a more consistent approach should be close to the one in 
the WS approximation in the case of the regular situation for the neutron spectrum. 

As it is seen in Table 4, at intermediate densities, fcp=0.6 -v- 0.8 fm~^, the regular situation 
takes place with the approximate equality Ap ~ Ainf for both kinds of boundary conditions and 
Ap'^^ ^ AP*-^^. Therefore one can expect that the gap function A„(r) will be approximately the 
same in the BCl and BC2 cases. At A;p=0.8 fm~^, where the equilibrium Z value is the same 
in both versions, these functions are displayed in Fig. 8. In this case the difference between 
Ap'-^^ and A^^^ is about 10%. As it can be seen in Fig. 8, the difference A^'^^(r) - Aj^*-^^(r) is 
also about 10% with the exception of small r < 3 fm which doesn't contribute appreciably to 
the matrix elements of A„. Thus, the accuracy of predictions for the gap function within the 
WS approach in the P2 model for the density interval under consideration is about 10%. Fig. 9 
collects the predictions for the gap functions for these values of kp which are taken in accordance 
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Figure 6: The neutron gap function A„(r) for kp=l.l fm~^, Z=20 and Z=24, in the BCl case 
(sohd hues) and in the BC2 one (dots), within the P2 modeL 
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Figure 7: The neutron single-particle spectrum for /cp^l.l fm ^, Z—20, in the BCl case 
(left) and the BC2 one (right), for the P2 model. 
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Figure 8: The neutron gap function A„(r) for A;f=0.8 fm ^, Z—AO, for the P2 model in the 
BCl case (sohd hne) and in the BC2 one (dashed hne). 



with the prescription of Section 3, i.e. those for the BCl kind of boundary conditions. 

Let us go to higher density values, kp > 0.9 fm~^, where the difference between the self- 
consistent values of Ap*"^ and Ap*"^ is significant. As the above discussion for /cf=1.1 fm~^ 
showed, there are two possibilities for the choice of A„(r) in this case with close results. To be 
definite, let us use the BCl version as the basic one if the relation A^^^ ~ Ajnf takes place. ^ 
In the opposite case, Ap*-^^ << Ajnf, we use the equilibrium Z^'^^ value and the gap function 
A^^^{r) for this Z. Corresponding gap functions are displayed in Fig. 10. The analysis similar 
to that for A;p=l.l fm^^ shows that the accuracy of the predictions for A„(r) is again about 
10% and only for A;p=1.2 fm~^ it is a little worse. 



^ There is a tiny difference between the values of Ajnf in the two versions of the boundary conditions which 
originates from a smaU difference of the two values of the asymptotic density. However, it is much less than the 
effects under discussion. 
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Figure 9: The neutron gap A„(r) for kp=0.6 ~- 0.8 fm ^ in the case of the P2 model. 

5 The P3 model 

The P3 model is similar to the P2 one, but now the many-body suppression factor in Eq. (0) 
is equal to /ni-b=l/3. The calculation scheme and the presentation of results is quite similar 
to that in the previous section. Values of the binding energy per a nucleon Eb in the P3 model 
are given in Fig. 11 for different densities, kp=0.6 -^ 1.2 fm~^, similarly to Fig. 5. Again the 
detailed analysis is made for kp=0.9 fm~^ and partially for /cf=0.7 fm~^. Comparison with 
Fig. 5 shows that positions of the absolute minima of the function Eb{Z) are the same as in 
the P2 model and in the P3 one for both kinds of boundary conditions at all kp values under 
consideration with only one exception, kp=0.9 fm~^ and the BC2 version. In the latter case 
the equilibrium Z values differ by 2 units. The main ground state characteristics of the neutron 
star inner crust within the P3 model are presented in Table 5 which is similar to Table 3 in 
Section 4 or Table 1 in Section 3. As one can see, the differences between all the values in Table 
5 and those in Table 3 are quite small. 

Let us now turn to the analysis of the neutron gap in the P3 model. The corresponding 
gap characteristics are presented in Table 6 which is similar to Table 4 (the P2 model) or Table 
2 (the PI model). Comparison with Table 4 shows that again the main effect is the general 
decrease of all the gap characteristics by the factor of 2/3 which corresponds to the values of 
the many-body suppression factor /m-b=l/3 in Eq. ((T)) in the P3 model and /m_b=l/2 in the 
P2 model. In two previous sections we discussed a pseudo effect of strong suppression of the 
neutron gap in some cases with high values of kp > 1 fm~^. It originates due to non- regular 



20 




5 10 15 
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Figure 10: The neutron gap A„(r) for A;f=0.9 -=-1.2 fm~^ in the case of the P2 model. 

behavior at big k-p of the neutron single particle spectra (the "pseudo Shell effect") in the 
WS approach. In the PI model (Section 3) the first case of a moderate suppression occurs at 
/cf=1-0 fm~^, in the P2 model (Section 4), at A;f=0.9 fm~^. In addition, in the P2 model, at 
k-p > 1.0 fm~^, the gap almost vanishes in some cases for the BCl case or the BC2 one. Table 
6 shows that in the P3 model this pseudo effect becomes even stronger, namely, the first case 
of vanishing occurs at A;f=0.9 fm~^. A typical case of such vanishing is shown in Fig. 12. The 
way to improve this drawback of the WS method and to find the neutron gap function A„(r) 
in such "bad" cases in the P3 model is the same as above. The predictions for A„(r) within the 
P3 model are displayed in Fig. 13 for /cf=0.6 -j- 0.8 fm~^ and in Fig. 14 for /cf=0.6 -j- 1.2 fm~^. 
The method to choose the version (BCl or BC2) at every value of kp is the same as it was 
suggested in the previous section for the P2 model. 
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Figure 11: Binding energy per a nucleon for various kp for the P3 model in the case of the BCl 
(sohd circles) and the BC2 (open circles) kinds of the boundary conditions. 
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Figure 12: The neutron gap for A:f=1-0 fm ^, Z=20 and Z=24, in the BCl case (sohd hues) 
and in the BC2 one (dots), for the P3 modeL 
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Table 5: Comparison of characteristics of equilibrium configurations of the WS cell for two 
different kinds of the boundary conditions in the case of the P3 model. 



kp, 
fm-i 


Z 


Rc, fm 


Eb, MeV 


/^n, MeV 


ID 1 




ID 1 


RP? 


Rn 

ID 1 


Rr'9 

ID vjZ 


0.6 


56 
54 


36.92 


36.96 
ou.uu 


2.3132 

Z.oloo 


2.3130 


3.5474 


3.5634 
o.ooio 


n 7 
u. / 


46 
48 


30.27 
31.09 


30.27 
30.85 


2.9603 
2.9630 


2.9658 
2.9636 


4.4009 
4.3857 


4.2937 
4.2913 


0.8 


40 


25.89 


26.27 


3.6016 


3.6079 


5.1069 


5.1355 


0.9 


20 
24 


18.40 
20.11 


18.62 
19.77 


4.2199 
4.2592 


4.2653 
4.2462 


6.5184 
5.9920 


5.8724 
5.6863 


1.0 


20 

24 


16.56 
18.25 


16.97 
17.85 


4.8642 
4.9576 


4.9581 
4.9001 


7.8971 
7.0551 


6.8332 
6.0549 


1.1 


20 
24 


15.05 
16.50 


15.45 
16.19 


5.5735 
5.7496 


5.7477 
5.6340 


6.2795 
8.3960 


7.9492 
6.7278 


1.2 


20 
26 


15.05 
15.30 


14.23 
14.90 


5.5734 
6.9276 


7.2771 
6.4566 


6.2792 
8.9663 


9.0045 
8.7410 



Table 6: Average gap characteristics in the P3 model 



kp, 
fm^^ 


Z 


kf, fm~^ 


A(0), MeV 


Aa., MeV 


Af, MeV 


Ainf, MeV 


MeV 


BCl 


BC2 


BCl 


BC2 


BCl 


BC2 


BCl 


BC2 


BCl 


BC2 


0.6 


56 
54 


0.5817 
0.5787 


0.5788 
0.5792 


0.721 
0.755 


0.755 
0.701 


0.719 
0.732 


0.713 
0.704 


0.723 
0.722 


0.715 
0.713 


0.778 
0.774 


0.774 
0.775 


0.81 


0.7 


46 
48 


0.6752 
0.6826 


0.6782 
0.6739 


0.789 
0.924 


0.850 
0.765 


0.832 
0.864 


0.852 
0.797 


0.806 
0.864 


0.865 
0.817 


0.893 
0.901 


0.896 
0.891 


0.92 


0.8 


40 


0.7676 


0.7791 


0.699 


1.019 


0.787 


0.923 


0.799 


0.943 


0.958 


0.964 


0.91 


0.9 


20 
24 


0.8489 
0.8747 
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Figure 13: The neutron gap A„(r) for kp=0.6 -r 0.8 fm ^ in the case of the P3 model. 

6 Discussion and Conclusions 

Recently a semi-microscopic self-consistent quantum approach was developed |2] for describ- 
ing the inner crust structure of neutron stars within the WS method with taking into account 
pairing correlation effects. It is based on the generalized energy functional method which 
is a modified version of the original Kohn-Sham one j3] for the case of superfiuid systems. In 
this approach, the energy functional is constructed by matching the realistic phenomenological 
functional by Fayans et al. j4j for describing the nuclear-type cluster in the center of the WS 
cell to the one calculated microscopically for neutron matter. The anomalous part of the latter 
was calculated within the BCS approximation. In this paper we take into account, in an ap- 
proximate way, corrections to the BCS theory which are known from the many-body theory of 
pairing in neutron matter. 

Unfortunately, up to now there is no consistent many-body theory of pairing in neutron 
matter. However, there exists a conventional point of view [0] that the BCS gap value is 
suppressed due to various many-body theory corrections significantly, by a factor between 1/2 
and 1/3. In the method developed in [HEl, the set of the neutron matter gap values A„(/cf) 
at the Fermi surface for the interval of < fep < 1.35 fm~^ is the only input to the microscopic 
part of the superfiuid component of the energy functional. In fact, we limit ourselves with the 
interval of 0.6 fm~^ < kp < 1.2 fm~^ in which the neutron pairing effects are expected to be 
larger. It is worth to note that the values of kp > 1.1 fm~^ should be considered as optional 
as far as the WS configuration with spherical symmetry is evidently unstable in this density 
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Figure 14: The neutron gap A„(r) for /cf=0.9 -t- 1.2 fm ^ in the case of the P3 model 
region [13113. 

We use a simple model to take into account approximately the many-body corrections to 
the BCS theory. In this model, the BCS value A^^^(A;f) is suppressed by a density independent 
factor which was taken to be /mb=l/2 in the first version of the model (named P2) and /mb=l/3 
in the second one, P3. These corrections influence the equilibrium configurations {Z, Rc) at 
different kp. The maximal variation from the BCS theory (the PI model) to the P2 version 
occurs at kp=0.7 fm~^, the equilibrium Z value changing by 6 units. As to the difference of the 
{Z, Rc) configurations found within the P2 and P3 models for the same version of the boundary 
conditions, it is usually negligible. The most important variation for P2 versus P3 occurs in 
the neutron gap function itself. We think that the realistic situation takes place somewhere 
between the P2 and P3 models. 
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